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ABSTRACT 

Using cosmological N-body simulations, we investigate the influence of the matter 
density parameter Q m and the linear theory power spectrum P(k) on statistical proper- 
ties of the dark matter halo population — the mass function n(M), two-point correlation 
function £(r), and pairwise velocity statistics vuir) and au(r). For fixed linear theory 
P(k), the effect of changing Vt m is simple: the halo mass scale M* shifts in proportion 
to f2 m , pairwise velocities (at fixed M/M*) are proportional to fij^ 6 , anc ^ na ^° cluster- 
ing at fixed M/M* is unchanged. If one simultaneously changes the power spectrum 
amplitude os to maintain the "cluster normalization" condition cigf^ 5 =const, then 
n(M) stays approximately constant near M ~ 5 x 1O 14 /i -1 M0, and halo clustering and 
pairwise velocities are similar at fixed M. However, the shape of n(M) changes, with 
a decrease of O m from 0.3 to 0.2 producing a ~ 30% drop in the number of low mass 
halos. One can preserve the shape of n(M) over a large dynamic range by changing 
the spectral tilt n s or shape parameter T, but the required changes are substantial - 
e.g., masking a decrease of Vt m from 0.3 to 0.2 requires An s ~ 0.3 or Ar ~ 0.15. These 
changes to P(k) significantly alter the halo clustering and halo velocities. The sensi- 
tivity of the dark halo population to cosmological model parameters has encouraging 
implications for efforts to constrain cosmology and galaxy bias with observed galaxy 
clustering, since the predicted changes in the halo population cannot easily be masked 
by altering the way that galaxies occupy halos. A shift in Q m alone would be detected 
by any dynamically sensitive clustering statistic; a cluster normalized change to as and 
VL m would require a change in galaxy occupation as a function of M/M*, which would 
alter galaxy clustering; and a simultaneous change to P(k) that preserves the halo mass 
function would change the clustering of the halos themselves. 

Subject headings: cosmology:theory — dark matter — galaxies: formation — large-scale 
structure of universe 
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1. Introduction 

In cosmological models dominated by cold dark matter (CDM) and an unclustered energy com- 
ponent (such as a cosmological constant), gravitational instability of primordial density fluctuations 
produces a population of dark matter halos, each in approximate virial equilibrium. Depending on 
their mass, these halos may host individual galaxies, galaxy groups, or rich galaxy clusters. Pro- 
vided that one focuses on systems of overdensity pip ~ 200, and regards higher density structures 
as subsystems within their parent halos, the halo population itself is insensitive to the gas pressure 
forces that influence the sub-dominant baryon component. In this paper, we investigate whether 
changes to cosmological parameters — specifically the matter density parameter Q m and parame- 
ters that describe the shape and amplitude of the primordial power spectrum — always produce 
measurable changes in the dark halo population, or whether two models with different combinations 
of these parameters can give rise to halo populations that are effectively indistinguishable. 

Our interest in this question is spurred by recent developments in the theory of biased galaxy 
formation. The uncertain relation between galaxies and dark matter is the primary limitation 
in testing cosmological models against observations of galaxy clustering. The "halo occupation 
distribution" (HOD) characterizes this relation statistically in terms of the probability distribution 
P(N\M) that a halo of virial mass M contains N galaxies of a specified type, together with 
prescriptions that specify the relative spatial and velocity distributions of galaxies and dark matter 
within these halos. Numerous recent papers have shown that the HOD framework is a powerful 
tool for analytic and numerical calculations of clustering statistics, for modeling observed clustering, 
and for characterizing the results of semi-analytic or numerical studies of galaxy formation (e.g., 
Kauffmann, Nusser, & Steinmetz 1997; Jing, Mo, & Borner 1998; Kauffman et al. 1999; Benson et 
al. 2000; Ma & Fry 2000; Peacock & Smith 2000; Seljak 2000; Berlind & Weinberg 2002; Bullock 
et al. 2001; Marinoni & Hudson 2002; Scoccimarro et al. 2001; Yoshikawa et al. 2001; White, 
Hernquist, & Springel 2001). In particular, Peacock & Smith (2000), Marinoni & Hudson (2002), 
and Berlind & Weinberg (2002) have argued that the HOD can be determined empirically from 
observed galaxy clustering, given an assumed cosmological model that determines the mass function 
and spatial and velocity clustering of the dark halo population. Empirical determinations of the 
HOD can provide insight into the physics of galaxy formation, and they may sharpen the ability 
of large scale structure studies to test cosmological models, since a model with an incorrect dark 
halo population may be unable to match the data for any choice of HOD (see Berlind & Weinberg 
2002 for further discussion). 

Suppose that we find a combination of cosmology and HOD that reproduces all aspects of 
observed galaxy clustering. Can we infer that the cosmology and the derived HOD are both 
correct? If there is another cosmological model that predicts an indistinguishable halo population, 
then the answer is clearly no, since the combination of this alternative cosmology with the same 
HOD would predict identical galaxy clustering. These cosmological models would be degenerate 
in the sense that they could not be distinguished by galaxy clustering data without relying on a 
predictive theory of galaxy formation (which might yield different HODs for the two cosmologies) . 
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Cosmological models that predict distinguishable halo populations could still be degenerate in this 
sense, if changes to the HOD can mask the differences in the halo populations; we will speculate 
on this point in §4, but we reserve a quantitative examination of it to future work. 

We focus our analyses on CDM cosmological models with Gaussian initial density fluctuation 
fields. Motivated by cosmic microwave background measurements (Netterfield et al. 2002; Pryke 
et al. 2002), we restrict our attention to spatially flat models with a cosmological constant or 
VL m = 1, though we demonstrate in passing that galaxy clustering data at z = probably cannot 
distinguish a flat model from an open model with the same Q m . The parameters that define our 
cosmological models are: the matter density parameter S7 m ; the normalization ag of the power 
spectrum P(k), which is the r.m.s. fluctuation of the linear density field filtered with a top-hat 
filter of radius 8ft. -1 Mpc ; and the shape of P(k), which we characterize by the spectral index 
n s of the inflationary power spectrum and the shape parameter T of the matter transfer function 
(see, e.g., Efstathiou, Bond, & White 1992). Although existing data impose constraints on these 
parameters, individually or in combination, here we allow each to vary independently over a fairly 
broad range, so that we can isolate the physical effects of the matter density, the amplitude of mass 
fluctuations, and the shape of the power spectrum on the resulting halo population. We concentrate 
entirely on the halo populations at z = 0, since this is where galaxy clustering data will be good 
enough to allow empirical HOD determinations in the near future. 

Using iV-body simulations described in §2, we measure (in §3) mass functions, two-point cor- 
relation functions, mean pairwise radial velocities, and pairwise velocity dispersions of dark matter 
halos. First, we consider models that have the same initial power spectrum P(k) but different 
matter density parameter Q m . We then move to cluster normalized models in which ag is changed 
to compensate the change in Cl m , so that the amplitude of the halo mass function is kept approxi- 
mately fixed at a cluster scale. Finally, we investigate models in which both the amplitude and the 
shape of P(k) are changed in order to match the amplitude and slope of the halo mass function at 
the cluster scale. We review our results and briefly discuss their implications for cosmological tests 
in §4. Our results overlap previous numerical studies of the halo mass function and halo clustering 
(recent examples include Jing 1998; Governato et al. 1999; Colberg et al. 2000; Jenkins et al. 2001), 
but they differ in the examination of controlled parameter sequences rather than specific cosmolog- 
ical models. Also, while some of these studies have focused on the large scale correlations of cluster 
mass halos, we devote considerable attention to the lower masses and smaller spatial scales that 
are important for non-linear galaxy clustering. Throughout the paper, we adopt a Hubble constant 
Hq = WOh kins" 1 Mpc -1 = 70 kms" 1 Mpc -1 (Freedman et al. 2001). 



2. Numerical Methods and Tests 

In this work, statistics are mainly based on simulations run with a Particle-Mesh (PM) N- 
body code (Park 1990). We run PM simulations with different combinations of density and power 
spectrum parameters. Each simulation follows the evolution of 200 3 particles in a periodic cube 
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200/i _1 Mpc on a side, with a 400 3 mesh to compute the gravitational force starting at z = 19 
and advancing to z = in 40 equal steps of expansion factor a. The mass of each particle is 
2.78 x lO^m/i^M©. We adopt the parameterization of Efstathiou, Bond, & White (1992) for the 
fluctuation power spectrum P(k). We divide our simulations into three categories corresponding 
to parameter changes in the matter density and the normalization and the shape of the initial 
fluctuation power spectrum. For each category, four cosmological models are investigated with 
(n m ,n A ) =(0.2,0.8), (0.3,0.7), (0.4,0.6) and (1.0,0.0). We generally regard the (0.3,0.7) model as 
the central one for comparison. In the central model, the spectral index is n s = 1.0, the spectral 
shape parameter is T = 0.20, and the r.m.s. mass fluctuation within spheres of radius 8/i _1 Mpc is 
o 8 = 0.9, which is consistent with the observed abundance of clusters (Eke, Cole, & Frenk 1996). 
For each model, we generate four independent realizations, and we use the dispersion among these 
realizations to estimate the statistical error bars on our clustering measures associated with the 
finite simulation volume. The total volume is comparable to that expected for the Sloan Digital 
Sky Survey (York et al. 2000) out to its median redshift z me d ~ 0.1. We use a FoF algorithm (Davis 
et al. 1985) to identify the dark matter halos, and, unless explicitly stated otherwise, the linking 
length is 0.2 times the mean inter-particle separation. With this linking length, the FoF algorithm 
picks out structures of typical mass overdensity p/p ~ 200 (Davis et al. 1985). 

The advantage of a PM code is that it is relatively cheap to run multiple simulations with a 
large dynamic range in mass. The limitation of a PM code is its moderate force resolution, but 
this limitation is not too serious for our purposes, since we only identify halos and do not attempt 
to measure their internal structure. We perform two tests to make sure that the PM resolution is 
adequate for our purpose. First, we make a comparison between PM simulations and simulations 
run with the GADGET tree code (Springel, Yoshida, & White 2001). Second, we perform a self- 
similar scaling test for the PM code. 

GADGET is a publicly available iV-body code that calculates particle accelerations by the 
hierarchical tree method of Barnes & Hut (1986). Its advantages over the PM algorithm are its 
spatial and temporal adaptability. The gravitational softening length, e, can be set significantly 
smaller than the initial particle grid spacing to increase force resolution, though if e is too small 
2-body relaxation effects will begin to influence the halo population. Particles have individual 
timesteps, which can vary continuously to ensure accurate time integration given the adopted force 
resolution. GADGET has fully periodic boundary conditions for comoving integrations via the 
Ewald summation method (Hernquist, Bouchet, & Suto 1991). 

We have compared the results of two simulations with 100 3 particles in a lOO/i -1 Mpc box, one 
evolved with the PM code using a 200 3 mesh (thus having equal resolution to the larger simulations 
on which our results will be based), and one evolved with GADGET using e = 0.3/i _1 Mpc. The 
initial conditions of the two simulations were identical to each other and used the parameters of 
the central model described above. Slices through these simulations are plotted for comparison 
in Figure 1. The structures formed in the two simulations are very similar at the level of detail 
discernible in this plot. The halo mass functions of the two simulations, also shown in Figure 1, 
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trace each other nearly exactly. The PM simulation produces slightly (~ 10%) fewer halos below 
10 13 /i _1 Mq, which is to be expected given its lower force resolution; the diameter of an overdensity 
200 sphere of 1O 13 /i _1 M0 is 1.04/i _1 Mpc, or about two PM grid cells. When the mesh size on 
this PM simulation was increased to 300 3 , the resulting mass function was consistent with the 
GADGET simulation even on these small mass scales. 

In addition to the halo mass function, the quantities important to our analysis are the cor- 
relation functions for both the mass distributions and the halo populations, which are shown in 
Figure 2. As with the halo mass functions, the matter correlation functions for the PM and GAD- 
GET simulations match very well. They begin to deviate at separations smaller than ~ 0.5/i _1 Mpc, 
or under one mesh cell. The halo correlation functions produced by the two methods are also con- 
sistent. Larger simulations and multiple realizations of each model are efficiently produced with 
the PM code, which takes up to an order of magnitude less computing time at the force resolution 
adopted here. 

We next perform a self-similar scaling test for the PM code. Assuming an Einstein-de Sitter 
universe (fi m = 1), we use a pure power-law power spectrum P(k) oc k~ ns with n s = — 1 (and no 
transfer function modulation) for the initial density fluctuations of the scale-free model. A 400 3 
force mesh is used to follow the gravitational evolution of 200 3 particles. The normalization of 
the power spectrum is set so that the nonlinear mass M* at the final output corresponds to 400 
particles. We define M* by the condition <j(M*) = 6 C = 1.69, where cr(M*) is the r.m.s. linear 
mass fluctuation on mass scale M*. The r.m.s. fluctuation is a = 1.0 at about 1/26 of the box 
size. If we identify this scale as 8/i -1 Mpc , the implied size of the periodic simulation cube is then 
207.62/i _1 Mpc at z = 0. Earlier outputs can be identified with higher redshifts or with larger 
volumes at z = 0. In the language of the redshift description, the simulation is evolved for 40 
time steps from z = 19, with 4 outputs (at redshifts z =1.00, 0.60, 0.25, and 0.00), when the 
corresponding M* values are 50, 100, 200, and 400 particles, respectively. The following results 
are averaged over four independent realizations of the simulation, and the error bars represent the 
dispersion among the four realizations divided by a/3 to yield the error on the mean. 

Efstathiou et al. (1988) performed the first study on scale-free cosmological models using N- 
body simulations and showed that many quantities in scale-free models have well-defined scaling 
behaviors. In an Einstein-de Sitter scale-free model, the cumulative halo mass functions at redshifts 
z\ and Z2 have the scaling behavior 



The left panel of Figure 3 shows the cumulative halo mass function at different redshifts in our 
scale- free simulations (solid curves) . The dotted curves in Figure 3 are derived by scaling the mass 
function at z = according to the above scaling rule. The scaling properties shown in the figure 
are generally very good. The high mass end of the mass function drops systematically below the 
self-similar scaling as M* increases, a likely result of the finite simulation volume, which supresses 



N(M > M 2 ;z 2 ) = 




(1) 



-6- 



the amplitude of fluctuations on the largest scales. Comparison of the first and last outputs suggests 
that this effect becomes noticeable above M ~ 2 x 10 15 /i _1 Mq at z = (for our adopted a$ = 1.0 
scaling), or a mass scale ~ 0.1% of the total simulation mass (10% in length scale). 

The iV-point correlation functions £(r) are expected to have a similarity transformation with 
the variable s = r(l + 2;) 2 /( 3 +™ s ) for an Einstein-de Sitter universe (Efstathiou et al. 1988). The 
scaling rule for the correlation functions at redshifts z\ and zi is then 

/l + Z2 \ 2 /(3+n s ) 

t{r 2 ;z 2 ) = n = I - - - - J r 2 . (2) 

From Figure 3 and Figure 4, it is clear that mass correlation functions and halo correlation functions 
at different values of M/M* obey the scaling rule quite well. There are some noticeable departures 
from self-similar scaling in the M*/2 panel at the two earliest output times (especially the earliest), 
when M*/2 corresponds to 25 and 50 particles, respectively. 

The clustering of halos is biased relative to that of the matter. The halo bias factor bh can 
be defined via ^(r) = b^ m (r), where £j, and £ m are the two-point correlation functions of halos 
and matter, respectively. The relation between the clustering of halos and that of the matter has 
been extensively studied based on analytical models and numerical simulations (e.g., Cole & Kaiser 
1989; Kashlinsky 1991; Mo k White 1996; Jing 1998, 1999; Porciani, Catelan, & Lacey 1999; Sheth 
& Lemson 1999; Sheth & Tormen 1999; Sheth, Mo, & Tormen 2001; and references therein). Mo 
& White (1996) give an analytic formula for the bias factor of halos of a given mass based on 
an extended Press-Schechter analysis (Press & Schechter 1974). This formula is quite accurate for 
halos with M > M*. Jing (1998) empirically modifies the formula to also fit the iV-body simulation 
results for halos of lower mass, 

, , \ (0.06-0.02n s ) , 9 _ i \ 

b h (M,z)={— A +l\ + u = 5 c (z)/a(M), (3) 

where a{M) is the r.m.s. fluctuation of the mass density at a mass scale M, 8 c (z) is the threshold 
density contrast for collapse of a homogeneous spherical overdense region at redshift z, and 5 Cj o ~ 
1.69 for an Einstein-de Sitter universe. For n s = —1 here, v = (M/M*) 2 / 3 . We derive the square of 
the bias factor numerically as a function of M/M* at different redshifts by averaging £h(r)/£ m (r) 
over comoving separations between 5/i _1 Mpc and 30/i _1 Mpc . In this range, the ratio is almost 
constant. We compare the computed bias factor with that given by Jing's fitting formula in Figure 5. 
Two conclusions can be drawn immediately from the comparison. First, since bias factors at 
different redshifts overlap with each other when plotted as a function of M/Af#, the simulation 
displays the correct scaling behavior. Second, the bias factors agree well with Jing's fitting formula. 
At higher mass, the computed bias is somewhat lower than that given by Jing's formula. In fact, 
the same trend can be found in Jing's plot (see his Figure 2). We attribute this slight deviation in 
part to a sample volume effect and in part to systematic errors in the approximate analytic formula 
given by Mo & White (1996) (see Sheth, Mo, k Tormen 2001). 
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Both the GADGET simulation test and the scale-free model test assure that the numerical 
artifacts in the PM simulations have almost negligible effect on the statistics we care about in this 
paper. For example, the mass scales that we are concerned with are generally higher than the 
level where the PM and GADGET halo mass functions begin to deviate slightly. Moreover, the 
correlation functions are measured at scales > 0.5 mesh cells, where the PM and GADGET results 
are consistent. The PM simulations reproduce the expected scaling behaviors in the scale-free 
model, which also demonstrates that we can measure masses of halos and their spatial distribution 
correctly based on these simulations. The scaling tests and code comparison suggest that there are 
some inaccuracies when the halo mass is lower than 50 particles or more than 0.1% of the total 
mass in the simulation, but even in these regimes there should be little impact on our conclusions 
because we compare different cosmologies evolved with the same code and numerical parameters. 
The PM code is therefore an ideal tool for our purpose in this paper. 



3. Comparison of Halo Populations 

3.1. Changing Q m with <jg, n s , and V Fixed 

We start with a simple case in which we fix the power spectrum (as linearly evolved to z = 0, 
with n s = 1.0, erg = 0.9, V = 0.20) and change Q rn . In linear theory, the evolved fluctuations 
in these models are identical, and this independence of VL m holds in the Zel'dovich approximation 
(Zel'dovich 1970) as well as its extension, the adhesion approximation, as emphasized by Weinberg 
& Gunn (1990). In fact, it holds remarkably well into the fully nonlinear regime, as shown in the 
top and middle panels of Figure 6, where slices of particle distributions at z = for models with 
different Q m but identical initial conditions are compared. These slices are taken from GADGET 
simulations that evolve the particle distribution from z = 19. Particle distributions of three spatially 
flat models with different values of Q m and one open model with Q m = 0.2 are extremely similar to 
each other (though particle velocities are higher for higher £l m ). Nusser k, Colberg (1998) extend 
the Zel'dovich/adhesion analytic explanation for this similarity of evolved structure by showing that 
the equations of motion of a collisionless gravitating system of particles in an expanding universe 
can be put in a form with almost no dependence on f2 m and if the linear perturbation growth 
factor is used as the time variable. The evolutionary histories of the four models in Figure 6 
are different, but this difference is captured almost entirely by the dependence of the linear growth 
factor on redshift. The bottom two panels demonstrate this point, showing the particle distributions 
of the flat Q, m = 1.0 and Q m = 0.2 models at the redshifts when the linear growth factor is 0.5, 
redshifts z = 1.0 and z = 1.74, respectively (we define the growth factor to be unity at z = 0.) In 
most comparisons of N-body simulations with different Q m , the amplitude and/or shape of P(k) is 
adjusted in concert with Q m . Such adjustments are well motivated by observational and theoretical 
considerations, but they mask the fact that the influence of Q m at fixed P(k) is extremely simple. 
Earlier examples of matched comparisons like those in Figure 6 include Figure 1 of Davis et al. 
(1985) and Figure 5 of Nusser & Colberg (1998). 
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Figure 6 suggests that halo mass functions for different models in this case should be nearly the 
same when measured as a function of particle number, i.e., mass divided by Cl m . This expectation 
is confirmed in Figure 7, which shows mass functions (from PM simulations) of four spatially flat 
models. In the left panel, all four (cumulative) mass functions have a similar shape and only shift 
horizontally relative to each other. When the halo masses are divided by O m , the four curves 
become nearly identical and overlap with the Q m = 1 curve. To allow a closer inspection, in the 
right panel we plot the ratio of the scaled mass functions to that calculated using Jenkins et al.'s 
(2001) fitting formula for the central model. The maximum relative difference between the four 
scaled mass functions is only about 10%. 

The 10% residual differences reflect slight changes in the characteristic densities and profiles 
of halos in different cosmologies. For Figure 7, we have identified halos with our standard FoF 
linking length parameter b = 0.2, which picks out structures with mean overdensity p/p ~ 200 
(Davis et al. 1985; Lacey & Cole 1993), close to the value p/p ~ 178 predicted for the post- 
virialization overdensity in the spherical collapse model for an Q m = 1 universe (Gunn & Gott 
1972; Peacock 1999). If we instead use b = 0.16, roughly doubling the overdensity threshold 
(0.2 3 /0.16 3 ~ 1.95), then halo masses drop by ~20% as previously linked particles are unlinked, but 
the level of agreement among different cosmological models stay the same. Another fairly common 
procedure is to scale b~ 3 in proportion to the virial overdensity predicted by the spherical collapse 
model, thus using a different b for each cosmological model. If we adopt this approach, taking 
the virial overdensities from Eke, Cole, & Frenk (1996), then the discrepancy of the scaled mass 
functions reverses sign and rises in amplitude to a maximum ~ 20%. We conclude that the effect of 
a pure fi m change is described by the simple scaling of halo masses to an accuracy ~ 10% for halos 
defined at fixed overdensity, that the discrepancies reflect the expected cosmology dependence of 
the characteristic virial overdensity, but that the magnitude of the cosmology correction predicted 
by the spherical collapse model is too large. Jenkins et al. (2001) also find better agreement among 
halo mass functions of different cosmological models for halos defined at constant overdensity rather 
than overdensities scaled according to the spherical collapse model. Henceforth we will use the fixed 
linking length b = 0.2 for all models, but our results would not be substantially different if we used 
b = 0.16 or scaled b with cosmology. 

The mass correlation functions £ m ( r ) for these four models are identical at large scales and 
gradually depart from each other at small scales (r < l/i -1 Mpc ), where the higher overdensities 
for lower £l m boost £ m (r). We calculate the two-point correlation function for halos of mass M by 
counting pairs of halos with masses in the range (M/ y/2, ^/2M) . Figure 8 shows halo correlation 
functions £h(r) for four values of M/M*. In each case, the four VL m models are indistinguishable 
from each other, as we would expect, since M* scales with £l m , and comparisons at fixed M/M* 
therefore remove any dependence on £l m . The difference in £ m (r) arises from differences in spatial 
structure within the halos themselves (which are difficult to discern in Figure 6). We also calculate 
the halo bias factor from the mass and halo correlation functions (not shown). When plotted as a 
function of M/M*, the bias factors of the four models are almost the same, as expected, and they 
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agree with Jing's (1998) fitting formula to a similar degree as seen in Figure 5. The halo correlation 
function plummets at small r because two halo centers cannot be separated by less than the sum of 
their virial radii. Halos for different Vl m have the same virial radius at a given M/M* (in simulation 
terms, the scaling of M* with Q m comes from the particle mass, not a change in halo size), so even 
this exclusion signature is the same in different models. 

Although mass functions and spatial clustering of halos are nearly identical as a function of 
M/M#, different values of O m produce different velocity fields. Figure 9 shows the mean pairwise 
radial (inward) velocity, vu = — (vi — V2) ■ (i"i — r2)/|ri — r2|, and the pairwise velocity disper- 
sion, a 12 = (v\ 2 ) — ( v i2) 2 , at M = M* and M = 8M*. Here Vj and r« (i = 1,2) are the velocities 
and positions of a halo pair, and the average is over all halo pairs with separations around r. Both 
v 12(f) an d 012(f) drop as tt m drops. The mean pairwise radial velocity 1712(f) increases for massive 
halos and increases much faster for pairs with small separations, a sign of the non-linear infall 
velocities induced by high mass halos. The velocity dispersion 012(f) continues rising out to large 
separations, and its amplitude at large scales seems to be nearly independent of halo mass. Linear 
perturbation theory predicts a relation between the peculiar velocity and matter density fields, 



where <5(x) is the mass density contrast and f(£l m ) ~ ( see > e.g., Peebles 1993). In the four 
cosmological models considered here, the mass density contrasts <5(x) should be the same, so linear 
velocity fields simply scale with f2^ 6 . The dotted curves in Figure 9 show the effect of dividing 
vi2(r) and (712(f) by f2^ 6 ; the agreement of these curves with each other and with the Q m = 1 curve 
demonstrates that the impact of a pure Q m change on halo peculiar velocities is well described by 
the linear theory scaling even on non-linear scales. 

Figure 6 suggests that the agreement of appropriately scaled halo mass functions, correlation 
functions, and pairwise velocity statistics should hold at other epochs for which the linear growth 
factors are equal. Though we do not show the results here, we have verified explicitly that the 
matching of scaled statistical measures in Figures 7, 8, and 9 holds equally well at the redshifts 
when the four models have growth factors of 0.5 (z = 1.0, 1.36, 1.50, and 1.74 for Q m = 1.0, 0.4, 
0.3, 0.2). Thus, for models with the same P(k) shape and the same as at z = 0, the influence 
of Q m and on the evolutionary history of the halo population is entirely encoded in the linear 
growth factor. 



Changing £l m on its own preserves the shape of the halo mass function but shifts the mass scale 
in proportion to VL m . This scaling implies that the threshold volume for collapse is independent of 
Q m , given the same power spectrum. If we require that the four cosmological models that differ 
in Q m have halo mass functions that match in amplitude at some physical mass scale (rather than 




(4) 



3.2. Cluster-normalized Models with Fixed P(k) Shape 
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M/M*), we need to adjust the power spectrum in a way that builds more massive halos in low-fi m 
models and fewer massive halos in high-f2 m models. In other words, we need a larger threshold 
volume for collapse in a low-f2 m model, so that its nonlinear mass M* increases. Therefore, if the 
shape of the power spectrum is kept unchanged, a higher ag is necessary in a lower Q m model in 
order to compensate for the lower mean mass density. 

Putting this idea in quantitative form, White, Efstathiou, & Frenk (1993) argued that the 
observed abundances of massive clusters of galaxies impose the constraint ag pa 0.57f2~°- 56 , with 
little dependence on the assumed shape of P(k). Many authors have revisited, refined and re- 
calibrated this "cluster normalization" constraint (see Pierpaoli, Scott, & White 2001 and refer- 
ences therein). One widely used formulation is that of Eke, Cole, & Frenk (1996), who found 
as = (0.52 ± 0.04)f2~ a52+ai3r2m for a spatially flat universe. Inspired by these results, we now 
consider a sequence of cluster normalized models in which we keep the shape of P(k) fixed but 
scale the amplitude as ag = 0.9(r2 m /0.3) -0 5 . 

Figure 10 shows cumulative halo mass functions for these cluster normalized models. Mass 
functions of models with Q m =0.2, 0.3, and 0.4 match in amplitude at a cluster scale (M ~ 
5 x 10 14 /i _1 M Q ). The mass function of the Q m = 1 model crosses the others at a somewhat lower 
mass scale, M ~ 2 x 10 14 /i _1 M Q . The discrepancy in matched mass scale just shows that the 
ag oc r^ ' 5 scaling does not hold all the way to Q, m = 1, and the Eke, Cole, & Frenk (1996) formula 
indeed implies that we should adopt a somewhat higher ag for Q m = 1. However, the crucial result 
for our purposes is that the mass functions of these cluster normalized models have systematically 
different shapes, and there is no value of ag we could choose that would make the Q m = 1 model 
match the others at more than a single mass scale. Over the smaller range Q m = 0.2 — 0.4, the 
cluster normalization condition keeps the halo mass functions reasonably well aligned over about a 
decade in mass, but by M ~ 1O 13 /i _1 M the halo space densities of the VL m = 0.2 and 0.4 models 
differ from those of the central model by ~ 30%, and from each other by nearly a factor of two. 

The reason that cluster normalization does not preserve the shape of the halo mass function 
is that raising ag increases M* but simultaneously drives down the space density of M* clusters 
by increasing the scale of non-linearity. For models with the same P(k) shape and the same value 
of M*, the space density of halos at fixed M/M* is proportional to VL m (see equation (7) below). 
The mass functions of cluster normalized models with different f2 m must therefore cross at different 
values of M/M*, and because the shapes remain the same as a function of M/M*, they cannot stay 
aligned as a function of M. 

In these cluster normalized models, the amplitude of £ m (r) increases as Q. m decreases, as shown 
in the upper left panel of Figure 11. At large scales, where the mass density field is still linear, 
£ m ( r ) is just proportional to the amplitude of its Fourier transform, the linear power spectrum, and 
hence to cr|. The roughly constant logarithmic offset in Figure 11 shows that this multiplicative 
scaling holds approximately even into the strongly non-linear regime. However, the remaining 
panels of Figure 11 show that the halo correlation functions are very similar at a fixed physical 
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mass. Qualitatively, this similarity makes sense, since a given mass scale corresponds to a higher 
value of M/M* when f2 m is higher, and the correspondingly higher halo bias tends to compensate 
for the lower mass clustering amplitude. However, it is impressive just how well this cancellation 
between mass clustering and halo bias works quantitatively. The halo exclusion signature does shift 
to larger scale for lower Q m because of the larger virial radii needed to compensate the lower mass 
density. 

Figure 12 shows that the mean pairwise velocity and pairwise dispersion increase with Q m in 
this cluster normalized sequence, but the dependence is weak, especially at large separations. As in 
Figure 9, more massive halos have higher pairwise velocities, especially at small separations. The 
weak dependence on Q m at large scales can be understood with the help of equation (4). In linear 
theory, velocity is proportional to f(fl m )as, which, when combined with the cluster normalization 
condition, reduces to approximately fi^ 1 . Non-linear evolution at small separations amplifies the 
difference between models, but only slightly. 

3.3. Matching Mass Functions 

By varying the amplitude of the power spectrum with Cl m , we can match the amplitudes of 
the cumulative halo mass functions at a given mass scale. However, this combination of parameter 
changes does not guarantee a match between the shapes of the mass functions at that mass scale. 
In order to match both the amplitude and the shape of the halo mass functions, we must adjust 
the relative amplitudes of fluctuations at different scales, or, in other words, the shape of P(k). 
We therefore need to determine the appropriate combinations of parameters (as and parameters 
related to the shape of the power spectrum) for each VL m model so that all models produce halo 
mass functions that agree both in amplitude and in shape. 

3.3.1. Analytic Solution to Parameter Combinations 

Instead of searching the parameter space to obtain the right power spectrum combinations, 
we make use of the analytic form of the halo mass function. There are three commonly used 
analytic formulae for halo mass functions. Press & Schechter (1974) developed the first theoretical 
framework based on the spherical collapse model (also see Bond et al. 1991) and provided the first 
theoretical formula (PS formula) for the halo mass function, which is still widely used. The PS 
formula is known to underestimate the halo abundance in the high mass tail (see Jenkins et al. 2001 
and references therein). Sheth & Tormen (1999) give an analytic formula (ST formula) obtained by 
fitting to the results of TV-body simulations, which turns out to be consistent with the ellipsoidal 
collapse model (Sheth, Mo, & Tormen 2001). Based on more simulation results, Jenkins et al. 
(2001) provide the most accurate fitting formula for halo mass functions. However, this formula 
includes an absolute value of a function, which makes it difficult to use for our present purpose: 
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obtaining an analytic solution of parameter combinations that lead to a good match in halo mass 
functions. We therefore use the slightly less accurate ST formula. 

In the notation of Jenkins et al. (2001), the mass function is defined in terms of the dimen- 
sionless function 

po diner 1 

and the ST formula for the halo mass function can be expressed as (Sheth & Tormen 1999) 

f{y) = A\f^- [1 + {av 2 yv} exp(-a^ 2 /2), v = S c /a(M), (6) 

V 7T 

where n(M) is the abundance of halos with mass less than M, S c ~ 1.69 is the threshold density 
contrast for collapse, a(M) is the r.m.s. fluctuation of the mass density at a mass scale M, po is 
the mean density of the current universe, A = 0.3222, a = 0.707, and p = 0.3. The mass function 
can be written as 

dn _ dn diner -1 _ p 
d\nM~ diner- 1 dlnM ""F 1 "'' [ ' 

where a = dlncr _1 /dlnM = (3 + n c g)/6 characterizes the local index n e ff of the power spectrum 
and is determined by the shape of the power spectrum. The slope of the mass function is 

dn \ dn / d In f d In a \ 



dlnM\dlnM/ dlnM \ d\av dlnM, 

We use a power-law to locally approximate the power spectrum (so that the last term of the 
right-hand side of equation (8) reduces to zero). 

We require that the halo mass function of a given cosmological model match as closely as 
possible that of the central model at some mass scale. In other words, at this mass scale, both the 
amplitudes (eq. 7) and the slopes (eq. 8) of the two halo mass functions should be equal. We thus 
solve these two equations for two unknowns (a,v) for each cosmological model. All the information 
contained in the power spectrum is fully embedded in (a,v). Therefore, by solving for these two 
parameters at a given mass scale, we can also determine the normalization as and shape of the 
power spectrum. 

The power spectrum can be expressed as P(k) oc k ns T 2 (k;T), where n s is the spectral index 
of the inflationary power spectrum and T(k; T) is the transfer function with shape parameter T. A 
change in either n s or V leads to a change in the shape of the power spectrum. We first keep V 
fixed and only tilt P (k). Since a = (3 + n e g)/6 and n e g is simply the sum of n s and an index given 
by the transfer function, it is straightforward to obtain the change in n s with respect to that of the 
central model: An s = 6Aa. Once the shape of P(k) is determined, the normalization er 8 can be 
obtained from v. Assuming that mass functions are matched against that of the central model at a 
cluster mass scale M = 5 x 10 14 /i _1 M©, we calculate n s and ag as a function of £l m . The left panel 
of Figure 13 shows that both n s and ag increase sharply towards low O m and drop slowly towards 
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high U m . In order to match with the halo mass function of the central model, we need n s =1.33, 
0.82, 0.42 and a 8 =1.16, 0.78, 0.55 for models with Q m =0.2, 0.4, and 1.0, respectively. 

Alternatively, we can keep the spectral index n s fixed (n s =l in our calculation) and change 
the shape parameter T. The right panel of Figure 13 shows T as a function of Q m when halo 
mass functions are matched at M = 5 x 10 h~ 1 M@. The shape parameter T also has a sharp 
increase towards low Q m and a slow drop towards high Q m . For models with Q. m =0.2, 0.4, and 
1.0, solutions to T are 0.35, 0.14, and 0.06, respectively. In this case, the normalization ag of each 
model is almost identical to that of the tilted model, which indicates that the two ways of changing 
the shape of the power spectrum are equivalent at this mass scale. Although pure-n s and pure-r 
changes alter the power spectrum in different ways, they produce power spectra (and thus mass 
functions) that match closely over a fairly large dynamic range, covering most of the regime in 
which a is of order unity. The near degeneracy can be broken by examining the low mass end of 
the halo mass functions. More generally, we could allow both n s and T to change. In that case, 
for each value of Q m , there would be a 1-dimensional locus in the n s — T plane along which models 
would satisfy the requirement of matching the central model's halo mass function at some mass 
scale. However, we expect that results from this more general case would be similar to those of the 
two simple cases. 



3.3.2. Simulation Results 

Using the power spectrum parameters given by the above analytic solutions, we run simulations 
and make comparisons with the central model. Here we only show results for different tilt models, 
since models with different values of V give similar results. 

Figure 14 shows cumulative halo mass functions for the four cosmological models. The mass 
functions agree with each other quite well over a large mass range (from ~ 10 13 /t _1 M Q to ~ 
1O 15 /i _1 M ), with only ~ 10% relative differences. We do not expect the cumulative mass functions 
to match exactly both in amplitude and in slope at M = 5 x 10 14 /i _1 M Q , for three reasons. 
First, the ST formula is not a perfect fit to the halo mass functions in the simulations. Second, 
we locally approximate the power spectrum as a power-law in order to obtain analytic solutions 
for the power spectrum parameters. Third, we require a match for differential mass functions 
instead of cumulative mass functions. Nevertheless, the cumulative mass functions trace each other 
remarkably well. The non-linear mass scales are quite different for these four models: M* = 
3.10 x lO^/T^M©, 1.03 x lO^/i^M©, 3.66 x lO^/i^M©, and 7.50 x W 8 h' 1 M Q , for Q m = 0.2, 
0.3, 0.4, 1.0, respectively. 

Figure 15 shows two-point correlation functions of the mass and halos. Mass is more strongly 
clustered in low Q m models, having a correlation function that is steeper and higher in amplitude 
than in higher £l m models. However, the clustering of halos exhibits the opposite trend. Halos in 
the Q, m = 1 model are highly biased with respect to those of other models, as one might guess 
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from the low nonlinear mass of this model, and this bias more than compensates for the weaker 
mass clustering. Conversely, at the mass scale M = 1.03 x 10 13 /i _1 Mq, shown in the upper right 
panel of Figure 15, halos in the Q m = 0.2 model are actually anti-biased with respect to the mass. 
Figure 11 shows that for cluster normalized models with fixed P(k) shape, the effects of bias almost 
exactly cancel the change in erg, leaving ^(r) nearly independent of Q m . Since our mass function 
matching approach imposes a redder P(k) shape (lower n s or V) for higher Q m , and increased large 
scale power amplifies the importance of bias, it is not surprising that bias over compensates mass 
clustering changes in these models, leaving £/i(r) steadily dependent on Q m . 

Figure 16 shows the mean pairwise radial velocity and pairwise velocity dispersion for the four 
cosmological models. As in the case of cluster normalized models, models with higher fi m have 
larger mean pairwise velocities and velocity dispersions than lower £l m models. The dependence 
on Q m is stronger and more systematic than that in Figure 12, again a sign of the redder power 
spectrum of higher Q m models. 



4. Summary and Discussion 

We have examined how changes to the matter density parameter Vt m and the shape and 
amplitude of the linear power spectrum P(k) affect the halo mass function, the two-point halo 
correlation function, and the first and second moments of the halo pairwise velocity distribution. A 
change in f2 m , at fixed P(k), simply shifts the halo mass scale. Therefore, if halo masses are scaled 
in proportion to Q m , halo populations of different Sl m models have identical mass functions and 
clustering properties. However, mean pairwise velocities and pairwise velocity dispersions, which 
scale as il^ 6 on large scales (and fairly far into the non-linear regime) break this degeneracy. A 
change of the power spectrum amplitude (<7g) along with Q m , maintaining a cluster normalization 
condition og oc f^°' 5 > produces halo populations that have nearly identical spatial clustering at 
a fixed physical mass scale and systematic but small differences in the mean pairwise velocities 
and pairwise velocity dispersions. However, while these cluster normalized models have the same 
halo space density at M ~ several xlO 1 M@, the shapes of their mass functions are systematically 
different. If we further allow the shape of the power spectrum to change (either by tilting it or by 
changing its shape parameter T), we can produce halo populations whose mass functions match very 
well both in amplitude and in shape, over a large mass range. However, the changes to the power 
spectrum shape cause correlation functions, mean pairwise velocities, and velocity dispersions for 
halos in the same mass range to differ. We conclude that the halo populations produced by distinct 
cosmological models are not degenerate. If they are indistinguishable by one statistic, they can be 
told apart using another statistic. 

Our results imply that a perfect knowledge of the halo population and its properties would allow 
us to pin down the underlying cosmological parameters - specifically the value of VL m and the shape 
and amplitude of P(k) - even before bringing in constraints from other cosmological measurements. 
However, a galaxy redshift survey detects galaxies rather than halos. Could changes to the halo 
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occupation distribution (HOD) mask the differences in the halo populations of different cosmological 
models? Models that differ only in £l m lead to the same halo mass function and halo clustering in 
terms of M/M*, but with M* oc f2 m . Their halos could be thus be populated the same way as a 
function of M/M* to produce indistinguishable galaxy spatial clustering. However, any dynamically 
sensitive statistics - virial masses of clusters and groups, pairwise velocity dispersions, the parameter 
combination f^ 6 jb inferred from redshift-space distortions, the galaxy-mass correlation function 
from galaxy-galaxy lensing - would distinguish models with different values of Q m immediately. 
Velocity bias within halos could mask some of these changes, but not all of them (Berlind k 
Weinberg 2002). Cluster-normalized models with different Q m and the same P (k) shape have 
different halo mass functions, and they would therefore require a different P(N\M) in order to 
keep the galaxy density fixed. This change would likely cause differences in the galaxy clustering 
even though the halo clustering is similar, since galaxy clustering is highly sensitive to P(N\M) 
(Berlind k Weinberg 2002 and references therein). Changing the shape of P(k) by the amount 
required to match mass function shapes makes substantial differences to the halo clustering and 
velocities, which seem unlikely to be masked by a change in the HOD. 

These speculations suggest that distinct cosmological models, which produce nondegenerate 
halo populations, cannot have galaxy populations that are indistinguishable in every spatial and 
velocity statistic, even if one allows complete freedom in the way that galaxies occupy these halos. 
We reserve a detailed investigation of this issue for future work. If our optimistic speculation 
holds true, then high precision measurements of galaxy clustering and galaxy-galaxy lensing should 
impose strong constrains on cosmological models without reliance on a priori models of galaxy bias. 

We thank Volker Springel for advice on GADGET and Changbom Park for the use of his PM 
code. DW thanks the Institute for Advanced Study for its hospitality and the Ambrose Monell 
Foundation for its financial support during the final stages of this work. This work was supported 
by NSF grant AST-0098584. 
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Fig. 1. — Comparison of PM and GADGET simulations evolved from the same initial conditions. 
Top panels show slices through the particle distributions, 100/i _1 Mpc on a side and 10/i -1 Mpc 
thick. Bottom panels show cumulative halo mass functions; in the right panel these are divided by 
the Jenkins et al. (2001) analytic mass function to allow more detailed comparison. The minimum 
group mass plotted corresponds to 10 simulation particles. 
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Fig. 2. — Matter and halo correlation functions of the PM (dotted lines) and GADGET (solid lines) 
simulations shown in Figure 1. The upper left panel shows the matter correlation functions and 
remaining panels show halo correlation functions for three different mass ranges, centered on M*/2, 
M*, and 2M*, where M* = 1.03 x 10 13 /i _1 M Q is the characteristic nonlinear mass (corresponding 
to 123 particles in these simulations). Results are averaged over four realizations, and error bars 
show the run-to-run dispersion divided by \/N — 1 = y/3 to yield the uncertainty in the mean. 
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Fig. 3. — Test of self-similar scaling in PM simulations with scale-free initial conditions. Solid lines 
in the left panel show halo mass functions at four evolutionary stages, when M*=50, 100, 200, 
and 400 particles (left to right). Dotted lines show the rightmost solid line scaled according to 
equation (1). Points in the right panel show two-point correlation functions of the mass at these 
four epochs. Dotted lines are derived by scaling the mass correlation function of the final output, 
when M* = 400M p , according to equation (2). Results are averaged over four simulations, and 
error bars on £(r) show the dispersion among four independent runs, divided by \/3- For clarity, 
error bars are plotted only for the stage of = 400M p , and error bars for other stages have 
similar magnitude. Physical scales are assigned by treating the final output as z = and adopting 
as = 1.0, which implies a comoving size L^ ox = 2O7.62/i~ 1 M of the simulation cube and a total 
mass M = 2.49 x 10 18 /i _1 M m in the cube. 



10 



AJtj? 0.1 



0.01 r 



a M,=400M 
- M. = 200M 
o M,= 100M 

• M. 50M p '. \ S< 
J i i i I i i 1 1 1 ll; 



-i — i — I i i i i I 1 1 — r 




1 5 10 

r (h _1 Mpc) 



5 10 

r (h _1 Mpc) 



50 



10 



: I I I I | 1 1 — I — | I I I I | 1 1 — r 



M = 2M, 



1 • 



Aj^ 0.1 



0.01 




I I I I I I 1_ 



10 F 



1 =- 



0.1 - 



0.01 - 



mm r 



i i i i I L 



1 1 1 I I I I | 1 1 rg 



M = 4M, 




_i i i i 1 1 1 1 



1 5 10 

r (h _1 Mpc) 



50 



1 5 10 

r (h _1 Mpc) 



50 



Fig. 4. — Test of self-similar scaling of the halo correlation functions, in mass ranges centered 
on M*/2, M*, and 4M*, as indicated. Points show halo correlation functions at the four 

evolutionary stages when M* = 50, 100, 200, and 400 particles. Dotted lines are obtained by 
scaling the result of the M* = 400 output (solid line) to earlier stages, according to equation (2). 
Error bars show the dispersion among four independent runs, divided by V3. 
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Fig. 5. — Halo bias factor as a function of M/M* in the scale-free simulations. Points show b 2 at 
the four evolutionary stages when M* = 50, 100, 200, and 400 particles. The solid line is calculated 
using Jing's (1998) formula. 
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Fig. 6. — Insensitivity of dark matter clustering to the values of Q m and ^a. The top and middle 
panels show slices of particle distributions at z = from simulations that have identical linear 
theory P(k) and Fourier phases but different combinations of £l m and as indicated. The two 
bottom panels are slices for models (1.0,0.0) and (0.2,0.8) at redshifts corresponding to a linear 
growth factor of 0.5. All slices have the size of lOO/i -1 Mpc x WOh' 1 Mpc x 10/tT 1 Mpc. All 
simulations were performed with GADGET using 100 3 particles in a (lOO/i" 1 Mpc) 3 volume and a 
force resolution of e = 0.3/i _1 Mpc. 




Fig. 7. — Influence of f2 m on the halo mass function. The left panel shows cumulative halo mass 
functions from simulations with four different combinations of (Q. m ,Q\) as indicated. A second 
curve for each low-f2 m model shows the result of dividing halo masses by £l m before computing 
N(> M). These scaled curves are almost perfectly superposed on the solid curve representing the 
Q m = 1.0 model, demonstrating that a change to Q m alone simply shifts the mass scale of the halo 
mass function. The right panel allows closer inspection of this result, suppressing dynamic range 
by dividing each mass function by the analytic fitting formula of Jenkins et al. (2001) computed 
for the parameters of the central model (fi m = 0.3, £l\ = 0.7). 




Fig. 8. — Influence of Q m on the halo correlation function, for halos in mass ranges centered on 
M*/2, M*, 4M*, and 8M* as indicated. Curves in each panel, nearly superposed, show £fe(r) for 
the four indicated combinations of (O m , ^a)- To preserve visual clarity, we plot error bars only on 
the central (0.3,0.7) model, but error bars for other models have similar magnitude. In each panel, 
open circles show the mass correlation function of the central model. 
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Fig. 9. — Influence of Vt m on the mean pairwise velocity (top) and the pairwise velocity dispersion 
(bottom) of halos with mass M ~ M* (left) and M ~ 8M* (right). Error bars are plotted only for 
the central model; they increase with VL m and are about twice as large for Q m = 1 as for Q m = 0.2. 
Dotted curves are obtained by dividing the velocity by S7^ 6 ; their agreement with the Q m = 1 
curve demonstrates that the Q m influence is captured almost entirely by the linear theory velocity 
scaling. 
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Fig. 10. — Halo mass functions for cluster normalized models, in which the shape of P(k) is fixed but 
the amplitude is scaled by oc J^ ' 5 - The left panel shows halo mass functions from simulations 
with four different (^l m ,Q\), as indicated. In the right panel, these mass functions are divided by 
the Jenkins et al. (2001) fitting formula for the central (0.3,0.7) model. 
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Fig. 11. — Mass and halo correlation functions for cluster normalized models. The upper left panel 
shows the mass correlation function. The remaining three panels show correlation functions of 
halos in mass ranges centered on M c , 4M C , and 8M C , where M c = 1.03 x 10 13 /i _1 Mq corresponds 
to M* of the central (0.3,0.7) model. Open circles show the mass correlation function for the central 
model. Error bars are plotted for the central model and have similar magnitude for other models. 
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Fig. 12. — Halo velocity statistics for cluster normalized models. Top panels show the mean pairwise 
radial velocities for halos with M ~ M c (left) and M ~ 8M C (right), where M c = 1.03 x 10 13 /T^Mq 
corresponds to of the central model. Bottom panels show pairwise radial velocity dispersions. 
Error bars are plotted for the central model and have similar magnitude for other models. 
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Fig. 13. — Changes in power spectrum parameters required to match the amplitude and slope of 
the halo mass function at M = 5 x 10 14 /i _1 M Q . The left panel shows the values of as (dotted 
line) and n s (solid line) required to match the mass function of a model with VL m = 0.3, ag = 0.9, 
n s = 1.0, and T = 0.2, when T is fixed at 0.2. The values of as we adopt for cluster-normalized 
models with fixed P(k) shape are plotted as a dashed line, for comparison. The right panel shows 
the required value of T if n s is held fixed at 1.0; the corresponding values of as are nearly identical 
to those in the left panel. Curves are calculated using the Sheth-Tormen (1999) analytic mass 
function, as discussed in the text. 
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Fig. 14. — Halo mass functions for models that have been matched at M = 5 x 10 h~ M 
by changing n s and as- As expected, the change in n s yields much better agreement of the mass 
functions than found for cluster normalized models with no change to the shape of P(k) (Figure 10). 
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Fig. 15. — Mass and halo correlation functions for the four models of Figure 14, in the same format 
as Figure 11. Halo mass ranges are centered on M c , 4M C , 8M C , where M c = 1.03 x l0 1? 'h~ l MQ 
corresponds to M* of the central model. Open circles in these three panels show the mass correlation 
function of this model. Error bars are plotted for the central model and have similar magnitude 
for other models. Changing n s to match the shapes of the halo mass functions leads to substantial 
differences in halo clustering, in contrast to Figure 11. 
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Fig. 16. — Halo velocity statistics for the four models shown in Figures 14 and 15, in the same format 
as Figure 12. Error bars are plotted for the central model; they increase with Q m and are about 
twice as large for £l m = 1 as for Vl m = 0.2. Changing n s to match the halo mass function shapes 
leads to significant systematic differences in the mean pairwise velocities and pairwise velocity 
dispersions as a function of Q m , larger than those in Figure 12. 
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